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SUMMARY 

General quenching-distance equations were derived as extensions of 
the original simple theory of quenching by diffusion of active particles. 
These equations take into account the possibility of gas-phase chain 
breaking and branching and the effect of the efficiency of the walls to 
destroy active particles. The general characteristics of the equations 
were examined. 

All available types of propane quenching data from the literature 
were treated in a consistent manner. The data covered a thirtyfold range 
of quenching distances and included the effects of temperature, pressure, 
and propane and oxygen concentrations. Lean-to-stoichlometric mixtures 
were treated. 

The effects of elevated initial temperatiire were predicted by means 
of the simple theory, and were found to agree qualitatively with the ob- 
served trends. 

The simple theory correlated all the available propane quenching 
data satisfactorily. From this agreement, it was concluded that gas- 
phase chain breaking and branching may be neglected in the treatment of 
quenching data for propane-oxygen-nitrogen flames. It was also concluded 
that the walls may be considered 100 percent efficient for the destruction 
of active particles. 

These conclusions are meaningfiil only in connection with the general 
concepts of flame quenching by diffusion of active particles, and do not 
establish the validity of the concepts . 


INTRODUCTION 

The quenching action of the walls on a flame is an important phe- 
nomenon in combustion. Consideration must be given to quenching in most 
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(3) Chain- branching reaction (The reaction of one active particle 
produces more than one new active particle.) 

(4) Chain-breaking reaction (The active particle is destroyed. ) 

Active particles diffuse to the walls of the quenching channel as 
well as into the unburned gas. It was assumed that a fraction e of the 
active particles that strike the walls are destroyed. The shape of the 
quenching channel has an effect on the amount of diffusion to the walls j 
it was assumed that the equations applicable to the case of a channel 
formed by infinite plane- parallel plates were approximately correct for 
long rectangular slots. Reference 8 investigated the end effects of rec- 
tangular slots and found that they were small in most cases. 

Straight- and branched- chain reactions promote flame propagation^ 
while chain breaking in the gas phase and at the wall tends to quench the 
flame. Thus^ the flame will be quenched if the chain- breaking reactions 
predominate. 


Derivation of Equations 

Reference 6 sets up the following simple criterion for a flame to 
propagate through a volume element in a quenching channel: 

ZviNi = AN^ (1) 

i 

Equation (l) states that a flame cannot propagate unless the total number 
of effective collisions of active particles is equal to a critical frac- 
tion A of the total number of molecules present in the gas phase. This 
criterion has been retained in the present work. Introduction of the 
factors of chain breaking and branching in the gas phase and of incomplete 
destruction of active par-tides at the walls merely changes the expression 
for the chain length v^. 

General equations; chain breaking and branching in gas phase; e/l . - 
In the derivation of the expression for the chain lengthy it is necessary 
to introduce chain-branching and -breaking coefficients f and g^ re- 
spectively. These coefficients have dimensions of reciprocal seconds and 
express the average frequency with which active particles take part in 
branching and breaking reactions. There is one solution for the case 
(f _ g) > 0 (branching predominates) and two solutions for the case 
(g - f ) > 0 (breaking predominates). The differential equation for dif- 
fusion of active particles in a duct formed by infinite plane-parallel 
plates is as follows (see refs. 12 and 13, and appendix): 
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d^Ci 


( 2 ) 


subject to the boundeiry conditions 

c . = c 




= ± d/2 


and. 


( 3 ) 


dc^ 


= 0 


dx 
X = 0 




The details of the derivation are given in the appendix. The following 
equations are obtained: 

Chain branching: 


d2 = 



( 4 ) 


where 

(f - g) >0, and d/2V(f - s)/^± < «/2 
Small amount of chain breaking: 


12 


d2 = 


A - 


(1 - e) 0.7 \ ^ 

l.KUl/ T 

kUf V 




kNf 

r Cl - e) 0.7 

V^i 

r pd2(g - f) 

/Toy 

^ k% 

/ ^i 

10 D° 

W J 


L ^ ■ g - f _ 

Z ! ^ 

i 




( 5 ) 


if MCA TN 5if09 


d 

width of rectangular slot^ cm 


diameter of cylinder^ cm 

f 

average chain-branching coefficient^ sec“^ 


chain-branching coefficient for active particles of one kind^ 

sec ^ ^ 

g 

N 

average chain-breaking coefficient^ sec”^ ^ 

Si 

chain-breaking coefficient for active particles of one kind^ 
sec”^ 

k 

average specific rate constant for reaction of 0^ and OH 

with fuel molecules^ cm^/(molecule)(sec) 

m 

power expressing variation of quenching distance with initial 
mixture temperature^ d oc Tq 

Nf 

number of fuel molecules per unit vol'ume^ number/cm^ 


number of fuel molecules per unit volume in unburned mixture 
at reaction temperature^ number/cm^ 

N. 

1 

number of active particles of one kind per unit volume^ 
number/ cm^ 


total number of molecules per unit volimie^ number/cm^ 

n 

power expressing temperature dependence of diffusion coeffi- 
cient^ D oc 

P 

total pressure of mixture^ atm 


equilibrium adiabatic partial pressure of 0^ and OH^ 

respectively^ atm 

Pi 

partial pressure of active particles of one kind^ atm 

Pi(eq) 

equilibrium adiabatic partial pressure of active particles 
of one kind^ atm 
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Tq initial mixture temperatvire, °K 

Tj^ mean reaction temperature in preflame zone^ °K 

TpCsi) equilibrium adiabatic flame temperatxire, °K 

x^y^z Cartesian coordinates 

6 fraction of total number of active particles striking walls 

that are destroyed 

chain lengthj that is^ average ntimber of effective collisions 
of an active particle of one kind with gas-phase molecviles 
before the particle is destroyed 



average time between effective collisions of an active particle 
of one kind with fuel molecules, sec 

f 

ACTIVE- PARTICLE DIFFUSION CONCEPT OF QUENCHING 
Description of Model 

The quenching equations were derived from a generalized form of the 
model of the quenching process described in reference 6. It was assumed 
that the light, rapidly diffusing active particles, H, 0, and OH, are 
necessary for flame propagation in hydrocarbon- air mlxtiores. The rate- 
controlling step in the combustion chain reaction is presiimed to be be- 
tween these active particles and fuel molecules and to take place ahead 
of the zone of maximum temperature rise. The concentration of active 
particles in the preflame zone is maintained by diffusion. This is the 
model treated by Tanford and Pease in their theory of burning velocity 
(refs . 10 and 11 ) . 

Foxir types of collisions were considered between active particles 
and other particles in the gas phase: 

(1) Collision without any essential effect upon the active particle 

(2) Reaction with a fuel molecule and subsequent regeneration after 

a chain of more rapid reactions (A series of such steps constitutes 
a simple, straight-chain reaction. ) 
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studies of other flame properties. For example^ Coward and Jones have 
pointed out the need to measure concentration limits of flammability in 
tubes large enough to preclude quenching effects (ref. l)j Lewis and 
von Elbe have described the relation of quenching by the electrodes to 
spark- ignition energy (ref. 2). 

For cylinders and plane-parallel plates^ the quenching effect may 
be characterized experimentally in terms of the quenching distance^ the 
smallest opening of a channel that will permit flame propagation. Quench- 
ing distance depends on the geometry of the channel and on the pressure^ 
temperature, and composition of the combustible mixture. 

Quenching offers a field for theoretical interpretation that has not 
been very thoroughly explored. In contrast to other flame phenomena, it 
is not difficult to assign plausible boundary conditions to the flame in 
proximity to walls. For example, the walls might be considered to be an 
essentially infinite sink for both heat and active particles (atoms and 
free radicals), the wall temperature to be the ambient temperature, and 
the active-particle concentration to be zero at the walls. 

References 2 and 3 describe a theory of quenching based on the loss 
of heat from the flame to the wall that utilizes the combust ion- wave 
theory of Lewis and von Elbe. The results were applied to quenching- 
distance data in references 2 and 4. Mixtures of methane and propane 
with oxygen and nitrogen were considered. A less comprehensive analysis, 
also based on heat loss to the walls, was presented in reference 5 to in- 
terpret the effects of various inerts on the quenching distances of 
hydrogen- oxygen- inert mixtures . 

More recently, equations have been derived that express the quench- 
ing distance in terms of a balance between the number of active particles 
lost by diffusion to the walls and the number that must diffuse into the 
fresh gas to sustain flame propagation (ref. 6). These equations have 
been tested with some success for a larger number of experimental situa- 
tions than either of the thermal theories. The equations in their pres- 
ent form apply only to lean-to- stoichiometric mixtures. Reference 6 
interpreted quenching distances for lean-to- stoichiometric propane-, 
ethene-, and isooctane-air mixtures at various reduced pressures} the 
geometrical factor to account for the difference between tubes and long 
rectangular slots was also predicted and reasonably well confirmed. Ref- 
erence 7 extended the application of the equations to quenching distances 
obtained for propane- oxygen- nitrogen mixtures over a range of pressures. 
Reference 8 dealt quite successfiolly with the predicted and observed geo- 
metrical differences between cylinders, annuli, and rectangular slots of 
various length-to-width ratios, as compared with the idealized case of 
infinite plane- parallel plates. 
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The equations developed in reference 6 were based on the following 
simple model of the chemical processes of flame propagation: (l) The 

rate- determining reaction is between fuel molecules and the light^ rap- 
idly diffusing active particles 0^ and OH; (2) there is no chain 
breaking or branching in the gas phase) (s) every active particle that 
strikes the walls is destroyed. 

The purpose of this report is to develop general equations for flame 
quenching^ based on the active-particle diffusion model of reference 6. 
The previous restrictions are removed. Chain breaking and branching in 
the gas phase are considered^ as well as incomplete destruction of active 
particles at the walls. The general equations are discussed in relation 
to all available types of quenching-distance data for lean-to- 
stoichiometric mixtures of propane with oxygen and nitrogen^ including 
data that show the effects of initial mixture temperature (ref. 9). All 
the data considered are for long^ rectangular quenching channels. The 
sensitivity of the results to the choice of diffusion coefficients and 
of the temperature dependence of the diffusion coefficients is also 
tested. 


<=1 
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SYMBOLS 

The following symbols are used in this report: 

fraction of total molecules present in the gas phase that 
must react if the flame is to continue to propagate 

concentration of active particles of one kind at any point 
across channel^ number/cm^ 

rate of production of active particles of one kind^ 
number/(cm^ ) (sec) 

concentration of active particles of one kind at the walls^ 
number/cm^ 

diffusion coefficient for active particles of one kind^ 
cm^/sec 

diffusion coefficient for active particles of one kind into 
unburned gas at atmospheric pressure and initial mixture 

temperature^ cm^/ sec 

diffusion coefficient for 0^ and OH^ respectively^ into 
unburned gas at atmospheric pressure and initial mixture 

temperature^ cm^/ sec 
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where 

(g - f ) > 0, and (g - f) d^/D^ < 1 
Large amount of chain breaking: 


d = 



( 6 ) 


where 

(g - f) > 0, and (g - f) > 2 

No chain branching or breaking in gas phase; e/l . - It was desired 
to derive the quenching equation for this less general case for direct 
comparison with the result derived in reference 6, in which gas-phase 
chain breaking and branching and incomplete destruction of active parti- 
cles at the walls were not considered. The differential equation of 
dif f us ion is : 


^2 

d 

®i 2 ~ ^ 

dx 


subject to the bo\ondary conditions 


and 


^i “ ^w, i 
X = ± d/2 




X = 0 


( 7 ) 


( 8 ) 


The details of the solution are given in the appendix, 
equation is obtained: 


The following 
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( 9 ) 


The solutions of equation (2) do not reduce to the solution of equation 
(j), when f - g is set equal to zero, because they are solutions of a 
basically different differential equation. 

It is also shown in the appendix that the geometry effect to be 
expected in comparing quenching in cylindrical and rectangular channels 
is the same as that predicted in reference 6, subject to the assumption 
£ = 1 . 


For the sake of brevity, following reference 7, equation (9) may be 
rewritten as 


where 



(1 - e) 0.7 \ ^ 
Ap L—J 
i 


Pi(eq) 



1 



Po 

Do 


poh\ 

Doh/ 


0.7 


( 10 ) 


( 11 ) 


APPLICATION OF QUENCHING EQUATIONS 
General Equations 

The general expressions for quenching distance (eqs. (4) to (6)) 
cannot be completely evaluated at present. The basic reason for this is 
that the mechanisms of active-particle reactions are not sufficiently 
understood. For example, the most general forms of the equations should 
allow for separate chain- branching and -breaking coefficients and 

g^ for each species of active particle. Furthermore, f^^ and gj^ are 

expected to vary with temperature, presstire, and mixtiore composition 
(ref. 12, p. 36). In view of the lack of information on these subjects, 
the specific coefficients have been replaced by average coefficients f 
and g in equations ( 4 ) to (6)j f and g are assumed to be the same 
for all active particles. A similar procedure was used in reference 6, 
where individual rate constants kj^ for the reaction between active 

particles and fuel molecules were replaced by an average rate constant k. 
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There is a similar lack of information with respect to t, the effi- 
ciency of the wall to destroy active particles. There have been many 
studies of the effects of wall materials on the low-temperatirre oxidation 
of hydrocarbons (ref. 2)^ but these results are probably not applicable 
to high-temperature oxidation in flames. Other studies specifically de- 
signed to measure the recombination of active particles on walls are also 
of doubtful value for the purposes of the quenching-distance equations. 

In such work^ care is taken to provide a known, reproducible sirrface in 
each experiment) in flame- quenching experiments, on the other hand, a 
wall that has been exposed to hot combustion products cannot be expected 
to remain uncontaminated. This is probably the reason for the usual 
failirre to observe any effect of different wall materials on quenching 
distance (ref. 9)j however, experiments in flame tubes, which were care- 
fiilly treated before the passage of each flame, did indicate that wall 
effects can be observed if the work is done properly (ref. 2, p. 335). 

In spite of the difficulties described, something may be gained by 
a discussion of the characteristics of the general quenching equations. 
Furthermore, ro\igh quantitative estimates of the quantities f - g and 
g - f may be derived from the equations, with the aid of calculated flame 
temperatures and partial pressures of H, 0, and OH. Such calculations 
have previously been reported (refs. 6 and 7) for the ranges of pressure, 
propane concentration, and oxygen concentration covered by the quenching- 
distance data of references 7 and 9. In order to include the data of 
reference 9 that show the effects of initial mixture temperature on 
quenching distance in the present discussion, some additional values of 
equilibrium adiabatic flame temperature and partial pressures of active 
particles were computed and are listed in table I. The calculations were 
carried out by the matrix method of reference 14, using the thermodynamic 
constants tabulated in reference 14 and the heat of formation of propane 
in reference 15. 

Inasmuch as the model contemplates a flame approaching the quenching 
channel and considers the events in some mean plane of reaction ahead of 
the flame front, it is necessary to choose values for the terms N^, Tp, 

and pj^. The choice made in reference 6 was to express the terms as 
functions of the equilibrium adiabatic flame temperature and active- 


particle concentrations, 
following relations were 

since these quantities 
chosen in reference 6^ 

can be calculated. The 
and are used herein: 


= 0.7 Tpi(eq) 

(12) 


Pi = 0.7 Pi(eq) 

(13) 


Nf = 0.5 %(Tj^) 

(14) 
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Chain branching; equation (4) . - A chain-branching reaction occurs 
if^ on the average^ each active particle that disappears in one of the 
- steps of the chain mechanism gives rise to more than one new active par- 

ticle. Thus, there is the possibility of a run-away reaction; in the 
limit, the chain length may become infinite, and cause a homogeneous ex- 
plosion. In such a case, the reaction cannot be quenched by the walls. 

The solution embodied in equation (4) is subject to the following 
^ restriction on the chain- branching parameter: 

CD 



Further, this is not merely a limit to the range of values for which the 
approximate expression for the chain length holds; it is also a physical 
limit. For the condition 



the chain length becomes infinite (see appendix). Since flames, not 
homogeneous explosions are of concern, and since the flames are at the 
quenching condition, the chain length clearly cannot be infinite. How- 
ever, equation (l6) may be used to estimate the upper limit for f - g. 

If is expressed in terms of the diffusion coefficient at 1 atmosphere 

and initial temperatiore. 



Substituting equation (l7) into equation (16) yields 


(17) 


- g) W - 1 

«2d0 \\J 


(18) 


Comparison of equation (l8) with the second bracketed term in the denomi- 
nator of equation (4) shows that the equation is consistent with the fact 
that homogeneous explosions are not quenched by walls ^ inasmuch as the 
condition for infinite chain length is also the condition for zero quench- 
ing distance. 


A mixture will certainly explode if the frequency of branching is 
the same as the frequency of react ionj that is^ if every reaction of an 
active particle produces more than one new particle. Now^ the average 
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rate constant for the reaction of active particles with fuel k seems 

to be -10” cubic centimeter per molecule per second (ref. 17 )j the 
frequency of effective collisions between active particles and fuel mole- 
cules kN|. is therefore 10^ to 10^ per second. If a single value of 

f - g is required to satisfy the criterion of equation (l8) for all the 
experimental points^ it should consequently be of the same order of mag- 
nitude as kN^ in the limit of infinite chain length. Inspection of the 

quenching- distance data of references 7 and 9 and of the calculated values 
of Tj^ shows that f - g may indeed be of the order of 10^ per second 

at the limit expressed by equation (l8). The limit of the equation is 
therefore consistent with the statement that a mixture will explode if 
the frequency of branching is the same as the frequency of reaction. 

It has already been stated that the coefficients f and g axe 
expected to vary with temperature^ pressure^ and composition of the mix- 
ture^ but that little or nothing is known about these variations in the 
case of hydrocarbon flame reactions. The quenching eqiiation (eq. (4)) 
can add nothing to the knowledge of the branching reaction itself. How- 
ever^ the term f - g^ expressing an over-all chain-branching reaction^ 
may be briefly discussed. If f - g is considered constant (despite the 
expected variations)^ the variables with which f - g is associated in 
equation (4) show that a constant amount of over-all chain branching 
would be a more important factor in the less vigorous mixtures in which 
large quenching distances are associated with relatively high pressures 
and low reaction temperatures. 

Chain breaking; equation (s) . - If a sufficient amount of chain 
breaking takes place in the gas phase^ no chemical reaction occurs in 
the limitj a flame cannot burn^ no matter how far the walls are separated. 

Equation (s) is one of the two solutions for the case of over-all 
chain breaking in the gas phase and is valid for small values of the 
chain-breaking parameter g - f . The restriction applicable to this 
solution is 


< 1 ( 19 ) 

■^i 

or, with written in terms of by means of equation (l7 ) ^ 

pd2(g - f ) /TqY 

D? W 


< 1 


( 20 ) 
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Equation (20) is not a physical limit; it is only a limit to the range of 
values of the chain-breaking parameter for which an expression for the 
chain length is a good approximation (see appendix). If g - f is re- 
quired to be a constant for all the propane quenching data^ its value is 

limited by equation (20) to about 10^ per second. This is, at most, only 
10 percent of kNf and is probably a negligible amount of chain breaking. 


It is seen from equation (s) that the quenching distance becomes in- 
finite under the following conditions: 


Pl^(g - f) ^ 1 

10 ■ 


( 21 ) 


Thus^ when this condition is satisfied^ gas-phase chain breaking is ca- 
pable of quenching the flame no matter how remote the walls may be^ and 
equation (s) is consistent with the general statement made at the begin- 
ning of this section. Equation (2l) is almost identical to equation (l8)^ 
which is the condition for zero quenching distance in the chain-branching 
case. As before^ the equation gives no information about the chain- 
breaking reaction itselfj but^ if a constant value of g - f is assumed^ 
it will be a relatively more important factor in weakly burning mixtures. 

Chain breaking; equation (S) . - The larger the fraction of molecules 
that must react to sustain flame propagation^ the smaller the amount of 
chain breaking that can be tolerated by the flame. Equation (6)^ which 
is applicable to larger values of the chain-breaking parameter g - f^ 
can be shown to take a form consistent with this statement. 

It will be shown in a later section that the assumption e = 1 is 
probably a good one. If equation (6) is rewritten subject to this assump- 
tion^ the following form is obtained: 



Here^ it is assumed as before that g - f is a constant. 
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Eqiiation (22) shovs that the quenching distance becomes infinite, 
that is, gas-phase chain breaking becomes overwhelming, if 


k% _ 

5 - fZ_y p 


(23) 


Since kNf> ~ 10^ to 10^ and ✓ ' Vj Iv ~ 10 "^^ then, very approximately, 

1 

(g - f) ~ i/a (24) 

Thus^ equation (6) leads to the expected relation between A and g - f . 


Treatment of Experimental Data 

All available types of quenching-distance data should be correlated 
by a suitable fonn of the equations^ if the diffusion concept of quenching 
is to be of value. The literature contains data that show the effects of 
pressure^ fuel concentration^ and proportion of oxygen in the oxidant 
atmosphere for propane- oxygen- nitrogen mixtures (ref. 7) and the effects 
of pressure, fuel concentration, and initial temperature, for propane-air 
mixtures (ref. 9). It should again be emphasized that the equations do 
not hold for rich mixtures in their present formj all the data to be 
treated in this report etre therefore limited to lean-to- stoichiometric 
mixtures . 

It has already been pointed out that the general equations (eqs. 

( 4 ) to (6)) contain the chain- branching and -breaking parameters f and 
g about which very little is known. Therefore, it is impossible to 
evaluate these equations numerically at present. It is, however, both 
possible and instructive to evaluate equation (lO), for the case of in- 
complete destruction of active particles at the walls, and the following 
equation, derived in reference 6 on the assumption that every particle 
striking the walls is destroyed: 


= ( 1 ^ i|r (25) 

The function ^ is defined by equation (ll). Equation (25) has been 
applied to propane quenching data previously (refs. 6 to 8)j but a con- 
sistent treatment has not been used for all the available types of data. 

Numerical evaluation of any of the forms of the equations involves a 
choice of the diffusion coefficients for the active particles into the 
unburned gas at Tq, and a choice of the temperature dependence of the 
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diffusion coefficients. This temperature dependence is expressed by the 
exponent n. Two different sets of values have been used in previous 
applications: (l) Reference 6 employed the diffusion coefficients cal- 
culated in reference 16, with a temperature dependence D (2) ref- 

erence 7 gave more sophisticated values of the coefficients, acco\mted 
for the small effects of propane concentration, and assigned a tempera- 
ture dependence D ^a,s decided in the present work to make 

separate comparisons of all the types of data for the two sets of diffu- 
fusion coefficients in order to determine whether the data were fitted 
better by one set of values than by the other. Table II summarizes the 
two sets of diffusion coefficients and the assigned temperature 
dependencies . 

The values of i|f calculated subject to a temperatiure dependence of 
the diffusion coefficients D oc T^ will be designated those cal- 

culated subject to Doc t 1*67 vill be designated \|f2. 

Limiting values of e; e^l (eg. (lO)) - - order to apply equation 
(10) to the experimental data, it was written in the following form: 



where 

X = (l - e) 0..7 (27) 


and 


= ^^yPH + Po + Poh) ^ 28 ) 

i 

The factor 0.7 (eq. (27)) expresses the average concentration of active 
particles in a volume element across the quenching channel ahead of the 
flame front and is introduced by eqiiation (l3). 

Several arbitrary values of the destruction-efficiency parameter K 
were chosen, and the quantity \|r (l - K2pj^(eq)/p) was calculated for both 

and •^ 2 ," bracketed term became negative in some cases if K was 

greater than about 10. Larger values of K were therefore riLLed out on 
physical groxinds, since they would lead to imaginary quenching distances. 
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This upper limit on K permits an estimate of the permissible range 
e. From equation (27), 


10 A 


0.7 


(29) 


The fraction of molecules A that must react for the flame to propagate 
refers to a plane ahead of the zone of maximum temperature rise on the 
fresh-gas side of the flame. It is therefore evident that A cannot be 
larger than the mole fraction of fuel at the lean limit of flammability. 

If A is taken as the lean limit ^ 0.0201 (ref. 17)^ equation (29) yields 
5 ~ 0.7. If^ on the other hand^ A is less than the lean limit^ e in- 
creases. Inasmuch as the upper limit is^ by definition^ 6=1^ the esti- 
mated range is 0.7 < e < 1. It should be noted that this result can only 
be considered approximate^ because it depends upon the choice of the fac- 
tor 0.7 to express the average concentration of active particles. 

Equation (26) does not correlate the quenching-distance data very 
well. A typical example is shown in figure 1^ a logarithmic plot of 

d^ against ^ 2 ^^ “ ^2pj^(eq)/p) for K = lOj it is assiomed^ as in refer- 

i 

ences 6 to 8^ that the ratio A/k is constant. It will be noted that 
there is a great deal of scatter in this correlation. 

Evaluation of simple quenching equations; 6=1 (eq. (25)). - Equation 
(25) is the quenching-distance equation developed in reference 6^ without 
consideration of chain breaking or branching in the gas phase^ and with 
the assiamption that every active particle striking the walls is destroyed. 
As stated previously^ equation (25) has been used to treat several types 
of quenching data (refs. 6 to 8)^ but the data of reference 9 that show 
the effects of initial mixture temperature have not previously been ana- 
lyzed. Therefore^ the agreement between the observed variation of quench- 
ing distance with initial temperature and the variation predicted by 
equation (25) is briefly described before proceeding to a general treat- 
ment of all the types of data. 

Values of ^ (eq. (ll)) were calculated for initial mixture tempera- 
tures of 0^^ 300^^ 400^^ 600^^ and 700^ K and a pressure of 0.9731 atmos- 
phere. The data of reference 9 were obtained for propane- air mixtures 
at the same pressure^ but for a more limited range of temperatures (300^ 
to 558^ K). The calculation of ijr for an initial mixture temperature 
of 0^ K was carried out with 300^ K as the base temperature for the dif- 
fusion coefficients in order to avoid making the denominator of equation 
( 11 ) equal to zero. 


Figure 2 shows plots of the observed quenching distance^ and 


and \|f 


1/2 


against initial temperature for four lean-to- stoichiometric 
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propane- air mixtures. in equation (25), the ratio 

to be inciependent of initial temperature^ the plots of 

against Tq show the predicted trend of quenching distance with initial 

temperature. Figiire 2 indicates that equation (25) gives qualitative 
agreement with experiment in regard to the order and relative spacing of 
the curves for 3.00^ 3.25^ 3.50^ and 4.03 (stoichiometric) percent by 
volume propane in air^ regardless of the assimiptions made for the diffu- 
sion coefficient. Furthermore^ the almost linear dependence of quenching 
distance on initial temperatiore observed over the relatively short range 
of conditions studied in reference 9 is predicted by the variation of 

and ^2 over the same range. 

The variation of quenching distance with temperat-ure may be expressed 
by the exponent m on the assumption that d oc Tq^. As reported in ref- 
erence 9^ m is about 0.5 for stoichiometric mixtures and increases for 
leaner mixtures. Values of m calculated by the method of least squares 
from the data of reference 9 are listed in table IIlj the values range 
from 0.60 for 4.03 percent by volume (stoichiometric) mixtures to 0.90 
for 3.00 percent mixt-ures. The predicted temperature variations may be 

1?2 -m 1/2 -m 

obtained from the analogous assiomptions^ oc Tq ^ and ^ *^0 

Values of m 2 ^ and m^ are also given in table and it is seen that 

the variations predicted from both and ^2 qualitatively correct 

with regard to the dependence on propane concentration. However^ the 
variation of *1/2 with Tq is less and that of tl/2 is more than is 

shown by the observed quenching distances. It should be noted that these 
observations concerning figiare 1 are governed by the assimption that the 
ratio A/k is independent of initial temperature; the results obtained 
when a reasonable temperature dependence is assigned are further discussed 
in a subsequent section. It may be concluded for the present that the 
simple diffusion concept of quenching gives fair predictions of the ef- 
fects of initial temperature. 

Figure 3 correlates all the available types of quenching-distance 
data by means of equation (25). Included are the effects of initial tem- 
perature^ pressure^ propane concentration^ and proportion of oxygen in 
the oxygen-nitrogen oxidant atmosphere (refs. 7 and 9). Other experiments 
that show only the effects of geometry (refs. 6 and 8) have been omitted. 

As stated previously^ only lean-to- stoichiometric mixtures are ^considered. 
Logarithmic plots are used in order to expand the scale; on straight plots^ 
many of the points would be near the origin. 

The left-hand portions of figures 3(a) and (b) show plots of d^ 
against and respectively. Straight lines with a slope of unity 


A/k is considered 
i|rl/2 and 
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are drawn through the points^ on the assijmption that the ratio A/k is 
a constant for all the data. It is seen that both and give 

good correlations j the included range of conditions causes about a thirty- 
fold variation in the quenching distance. 

However, the points for elevated initial temperatures deviate con- 
sistently in figure 3(a) from the line defined by the rest of the points. 
In figure 3(b), the corresponding points lie closer to the line. It 
therefore appears that increase of the initial mixtiire temperatiire above 
300° K introduces effects that are reasonably well accounted for by values 
of based on the diffusion coefficients of reference 7 with a tempera- 
ture dependence D However, close comparison of figures 3(a) 

and (b) shows that the points other than those for elevated temperatures 

are slightly better distributed about the line for d^ against \lr^. In 

addition, the possible dependence of A and k on initial temperature 
sho\ald be considered. 

Increase of the initial mixture temperature from 300° to 600° K, the 
approximate range of the data of reference 9, causes a maximum increase 
in the calculated mean reaction temperature T|^ of 200° K. Inasmuch as 

the data remaining correlate so well in figure 3, even though they rep- 
resent a range of from 1350° to 2000° K, the average rate constant 

is apparently insensitive to temperature variation fot the purposes of 
correlation. 

The constant A, defined as the fraction of molecules present in the 
gas phase that must react if the flame is to propagate, may on the other 
hand depend on the initial temperat\u:e. Reference 6 showed that quenching 
distances could be calciilated from equation (25) with an average deviation 
of only 3 percent when A was taken as the mole fraction of fuel at the 
lean flammability limit, and k was taken as the semiempirical rate con- 
stant derived from burning velocity measurements (ref. 16). If A is 
actually the mole fraction at the lean limit, there definitely should be 
a dependence on initial temperature) it is well known that the lean-limit 
concentration decreases nearly linearly with increased temperature (ref. 
l). Wohl has objected to such an interpretation of A (ref. 18 )j never- 
theless, A may be expected to vary with temperature in the same manner 
as the lean limit, even though it may not be correct to set it equal to 
the lean limit. This conclusion is drawn from the behavior of lean-limit 
flames: the increased enthalpy of a combustible mixture due to a higher 
initial temperature is just offset by an equivalent decrease in the heat 
of combustion of the limit mixture, so that a certain minimum enthalpy 
per unit area appears to be necessary for a flame to propagate (refs. 1 
and 19). It is reasonable that similar considerations should also apply 
to other than limit flames . 
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The factor A (eq. (25)) was therefore set equal to the mole fraction 
of fuel at the lean flammability limit. The value of A was taken as 
0.0201 for all propane- oxygen- nitrogen mixtures at 300° K (ref. 17). No 
data were foiind in the literature on the variation of the lean limit of 
propane with initial temperature. It was therefore assumed that the per- 
centage change in lean-limit equivalence ratio of propane with temperatirre 
is the same as that determined for n-pentane (ref. l). This assiimption 
leads to a predicted decrease of 0.19 percent by volme in the lean-limit 
concentration per 100° C temperature increase. The result is corroborated 
by a limited amount of unpublished NACA data^ which give the variation as 
0.20 percent by volume per 100° C temperature increase. 

Plots of d^ against A\lc 2 _ and A^lfg are shown in the right-hand 

portions of figures 3(a) and (b), respectively. It is again emphasized 
that variation of A with initial temperature is the important feat\ire^ 
not the absolute magnitude of A. The fig\rres show that the elevated- 
temperatiire data are correlated considerably better by Ailr^ than by 

the opposite is true for ilfg and A\lf 2 ^ because ijfg alone correlates 
these data quite well. 


DISCUSSION 

Figiare 3 shows that the simple diffusion theory of quenching corre- 
lates a wide range of quenching-distance data quite well. Variations in 
pressure^ temperature, propane concentration, and oxidant atmosphere are 
included. Similar resiilts are obtained with either of the two sets of 
diffusion coefficients. However, the logical assignment of a temperature 
variation for A somewhat favors a temperature dependence of the diffu- 
sion coefficients of the form D cx T^ rather than D oc jl-S? for pur- 
poses of the theory. 

It is suggested that quenching distances for propane- oxygen- nitrogen 
mixtures may be reasonably well calculated by means of equation (25) if 
A is taken as the lean flammability limit and k is taken as the value 

derived from the right-hand line of figure 3(a), namely, 2.25x10“^^ cubic 
centimeter per molecule per second. The value of A should be the ap- 
propriate one for the initial mixture conditions. 

The general quenching-distance equations developed in this report 
cannot be evaluated numerlcaLlly, because of insufficient knowledge about 
the chain- branching and -breaking parameters f and g. However, as 
previously pointed out, the equations show that, for constant values of 
f - g or g - f, these factors will be more important in the less vig- 
orous mixtures. Reference 7 Includes data for both weak mixtures 
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( 02/(02 + N 2 ) = 0,17) and vigorous mixtures ( 02/(02 + ^ 2 ) = 0.70)j inas- 
much as figure 3 does not show any consistent deviations of these data 
from the simple theory^ it is concluded that gas -phase branching and 
breaking need not be considered in the quenching of lean-to -stoichiometric 
propane -oxygen -nitrogen flames. 

Evaluation of the effects of incomplete destruction of active parti- 
cles at the walls (eq. (lO) and fig. l) shows that the efficiency of de- 
struction e cannot be very small. Comparison of figures 1 and 3 sup- 
ports this conclusion; the simple theory^ subject to the assiomption 
e = 1^ correlates the quenching data much better than does equation (lO). 

This discussion and the conclusions drawn are based on quenching data 
for propane- oxygen- nitrogen mixtiores. It is believed that the conclusions 
are valid for other hydrocarbon fuels as well^ inasmuch as propane con- 
tains both primary and secondary carbon-hydrogen bonds^ and is therefore 
representative of unbranched alkanes in general. However, it should be 
emphasized that the conclusions are meaningful only in connection with 
the general concepts of flame quenching by diffusion of active particles 
and that no proof has been advanced for these concepts. 


CONCLUSIONS 

General equations for quenching distance have been derived, based on 
the active-particle diffusion concept of quenching. The equations include 
the effects of gas-phase chain branching and breaking and of incomplete 
destruction of active particles at the walls. 

All the available types of propane- quenching data for long rectangu- 
lar slots may be well correlated by the simple theory, in which gas-phase 
chain breaking and branching are not included and complete destruction of 
active particles at the walls is assumed. The data represent about a 
thirtyfold variation in quenching distance, and include the effects of 
initial temperature, pressure, proportion of oxygen in the oxidant atmos- 
phere, and propane concentration; however, only lean-to- stoichiometric 
mixtiires are treated. 

Data on the effects of initial temperature on quenching distance, 
which have not previously been treated, agree qualitatively with predicted 
trends . 

Application of the simple theory to the data indicates that the tem- 
perature dependence assigned to the diffusion coefficients for active 
particles has some effect on the degree of correlation. For piirposes of 

the theory, a dependence greater than D cx is diffusion coef- 

ficient; T is temperatiire) is favored; however, this should not be taken 
as evidence for the true temperature dependence. 
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From the agreement of the data with the simple equation that ignores 
the complicating features^ it is concluded that chain branching and break- 
ing may be neglected in the treatment of quenching data for propane- 
oxygen-nitrogen systems. It is also concluded that the assumption of 
complete destruction of active particles on collision with the walls is 
valid. It should be emphasized^ however^ that the conclusions are mean- 
ingful only in connection with the general concepts of flame quenching 
by diffusion^ and that no proof has been advanced for these concepts. 

Lewis Flight Propiolsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland^ Ohio^ December 20^ 1954 
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APPENDIX - QUENCHING BY PLANE- PARALLEL PLATES OF INFINITE EXTENT 

For the case in which total destruction of active particles does not 
occur at the walls and in which chain breaking and branching may occirr in 
the voliome^ the diffusion equation describing quenching by plane-parallel 
plates of infinite extent is given by (ref. 13) 


7-Dj^Vcj^ + (f^ - g^) CjL + Co^i = 0 


(Al) 


where ^ for the i 
V • D . V c 4 


th 


species of active particles^ 


rate of disappearance of active particles (per unit 
volume) due to diffusion to walls 

Cq ^ rate of creation of active particles (per unit volume) 

^ by straight- chain reactions 

0 < (f^ - gi)cj_ rate of creation of active particles (per unit volume) 

due to branched- chain reactions 

0 < (g^ - rate of destruction of active particles (per unit 

volimie) due to chain- breaking reactions 

If it is assumed that average^ constant values may be used for 
equation (A1) may be written^ in rectangular coordinates^ 


/ d^c^ 

1 \~~z' ^ ~~T ) + (^1 " + '= 0 , 1 = 0 

\dx dy dz / 


(A2) 


Application to the case of the volume contained between plane- parallel 
plates of infinite extent and of separation d simplifies equation (A2) 
to 




subject to the boundary conditions 


(A3) 
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‘"i i 'N 


= ± d/2 


and. 




dc. 


= 0 


(A4) 


djx 

X = 0 ^ 

The three cases considered are 

(1) Chain-branching reaction, (fj^ - g^) >0 

(2) Chain-breaking reaction, (gj^ - fj^) > 0 

(3) Straight- chain reaction, - gj^) = 0 

Over-all chain branching, (f - gi) ^ 0« - Integration of equation 
(as) gives (ref. 13) 


Cj^ = A* cos 




X + B' sin 




'oj- 




(A5) 


Application of the boundary conditions (eq. (A4)) then yields 



- 1 + c 


w, i 


(A6) 


The active-particle chain length for a particular species is given by 


''i = 


'd/2 

f-d/2 


c^dx 


o, 1 1 







O, 1 1 


(A7) 
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Equation (A7) remains finite as long as ^ ^ For this range 

of values^ Vj^ can be closely approximated by (ref. 13) 


H 

C) 

1^ 

H* 

8d2(fj^ - g^) 

1 

+ 

pfl - Si)'’'! ^o, i^ij 

jt^Di 

- gi) 

S,i'"i 


From reference 6 


(A8) 


Pi 


- A V ""i 


(A9) 


Also 


Tj = 




(AlO) 


and 


D. = 

v^oy b. 


(All) 


Further, 


=0,1 = + (fi - 8i)u] 


(A12) 


and it is assumed that 


kj^ = k 


A 


= f \ for all i'i 


S± = & j 


and 


Cv^i - (l " e) 0.7 Ni 


(A13) 
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It should he noted that the form chosen for ^ represents the 
simplest one that satisfies the criteria c„ ^- = 0 for e = 1 and 
= 0.7 for e = 0. 

Using equations (AS) to (A13) yields 


d^ = 




(1 - e) 0.7 
. i 

Pi’ 

P 





" " k% 

+ (f - g 

) J 

\8 ho) 


k% 

r , (1 - 0 o.7i 

El 


1 


_L T" 

kNf 

1 _ 

pd2(f - g) 

(Toy 



^ f - g _ 

1 

_L — 

2 TiO 

rt 

Vr) 


(A14) 


Over-all chain breaking, (fi - g^) < 0. - Integration of eq^uation 
(A3) gives (_ref. 13) 


Cj^ = A" cosh _ — - X + B" sinh \ / - --^ — - x + — 


»i 


»i 


Si " 


(A15) 


Application of the boiondary conditions (eq. (A4)) then yields 


cosh 




Ca = 


ls± - fj^ 

®i 


+ C, 


cosh 


icr f ; 1 

/Si " ^i d ' ’ 


(A16) 


The active-particle chain length for a given species is given by 


v-i = 


‘d/2 

-d/2 


c^dx 


C_ 


' 0,11 


Thus , 


Vi = 


Ui - ^‘i)'^i 


‘^o, i'^i 


tanh 

d 

2 


— — ■ ■ I I ^ — .....z. 


ISi - fi 
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For values of (gj^ - fj^) equation (A17) may be very closely 

approximated by 


Vi = 


— 1 
H 

1 

O 

< 

H- 

ti^CSi - fi) 

" ^"Csi - fi)“ 


Ui - t^)\ Cc,i\ 

12 

10 

I 

<=o,i'^i 


(A18) 

If, on the other hand, (gj^ - fj^) d^/D^^ > 2, equation (A17) can be very 
closely approximated by 


''i = 


^w^l 


CSi “ ^o, 


a /Si - 


(A19) 


First, for the case where (gj^ - f^) d^/Dj^ < 1, equations (A9) to (A13) 
may be used, together with equation (A18), to show that 


12 




d^ = 


(1 - e) 0.7 


A - 


1 - 


g - f 
k% 


kM 


(1 - e) 0.7l 


pd^(g - f) 

^ kM 

^ - 1 


1 

H 

O 

g - f 

i 



(A20) 


Similarly, for the case where (gj^ - fj^) d^/Dj^ > 2, expression (A19) may 
be used to show that 


kU. 


f_ (l - e) 0.7 
1 


g - f ^ g - f 


kNf 


2Pi 


(g - f) P" 


w 


d = 




kNx 


(A21) 


/ / P \g - f/ 
i 


- A 
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Straight- chain reaction, - gj^) = 0. - In the case of a straight- 


chain reaction, equation (A3) reduces to 



(A22) 


Solution of equation (A22), subject to the boundary conditions (eq. (A4)), 
yields 


Setting and using equations (A9) to (All), (A13), and 

(A24), together with the assumption that = k for all i's, the fol- 
lowing equation is obtained for the quenching distance: 



It is interesting to note that the geometry effect, given in ref- 
erence 6 for plane- parallel plates and for cylindrical tubes, remains 
unchanged as a res\ilt of the assumption of nonzero active- particle con- 
centrations at the wall. This follows from the fact that solution of the 
diffusion equation for cylinders 




The active-particle chain length for a given species is given by 



(A24) 



(A26) 


subject to the boundary conditions 
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Cj^ = at the walls 


(dCj^/dr) = 0 at r = 0 
yields, for the chain length. 


d§ 


•=^1 32 


1 

^o, i'^i 


(A27) 


(A28j 


Thus, if is taken to he the same for both plane-parallel plates 

and cylindrical tubes, the geometry effect given in reference 6 remains 
unchanged . 
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TABLE I. - CALCULATED PARTIAL PRESSURES OF ACTIVE PARTICLES AND 
EQUILIBRIUM ADIABATIC FLAME TEMPERATURES 


Propane-air mixtures at 1 atm and various initial temperatures. 


Propane 
in air, 
percent 
by vol- 
ume 

Initial 
temper- 
ature , 
°K 

Partial 
pressure 
of H atom, 
atm 

Partial 
pressure 
of 0 atom, 
atm 

Partial 
pressure 
of OH radi- 
cal, 

atm 

Equi- 
libri- 
um flame 
tempera- 
ture, 

°K 

3.00 

0 

0. 01x10"^ 

0.16x10“^ 

2.77xl0"'^ 

1733 


300 

.10 

1.03 

9.64 

1941 


400 

.22 

1.76 

13.49 

2011 


600 

.85 

4.60 

25.17 

2148 


700 

1.53 

7.03 

33.10 

2214 

3.25 

0 

0.04x10"^ 

0.39X10"^ 

5.33X10"^ 

1843 


400 

.65 

3.06 

20.52 

2108 


600 

2.00 

6.89 

34.63 

2233 


700 

3.27 

9.85 

43.53 

2293 

3.50 

0 

0.14X10"'^ 

0.77X10-^ 

o.oexio""^ 

1949 


300 

.96 

3.06 

22.12 

2134 


400 

1.62 

4.46 

27.87 

2193 


600 

4.06 

8.88 

43.24 

2304 


700 

6.06 

12.09 

52.49 

2357 

4.03 

0 

1.50x10"^ 

1. 01x10"^ 

13.19X10“^ 

2122 


300 

4.91 

3.61 

28.90 

2273 


400 

6.37 

4.85 

34.20 

2313 


600 

11.55 

9.18 

50.11 

2401 


700 

15.19 

12.26 

59.52 

2444 

4.50 

0 

2 . 80 x 10 "'^ 

0.08X10"^ 

3.66x10*^ 

2073 


400 

10.82 

1.73 

20.95 

2312 


600 

18.60 

4.88 

37.92 

2414 


700 

23.79 

7.40 

48.19 

2462 

5.00 

0 

2.14x10"^ 

0.01x10“^ 

0.91x10“^ 

1966 


300 

7.34 

.12 

4.83 

2161 


400 

10.51 

.29 

7.87 

2229 


600 

20.40 

1.30 

18.84 

2358 


700 

27.36 

2.50 

27.32 

2419 

6.00 

0 

0.68x10"^ 


0.06X10'^ 

1762 


300 

3.18 


.50 

1958 


400 

5.07 

— 

.93 

2029 


600 

12.26 

0.05X10"^ 

3.05 

2175 


700 

18.47 

.14 

5.26 

2249 
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TABLE II. - DIFFUSION COEFFICIENTS 


FOR ACTIVE PARTICLES 


Diffusion coefficients 

a^l 

^2 

Dj£, cm^/sec (300° K, 1 atm) 

1.80 

1.02 

Dq, cm^/sec (300° 1 atm) 

.40 

.29 

Dq£, cm^/sec (300° K, 1 atm) 

.28 

.26 

n, (DcxrT^) 

2.00 

1.67 


®'Ref . 16 . 

^Diffusion coefficients into unburned 
stoichiometric (4.03 percent by volime) 
propane-air mixture} ref. 7. 


TABLE III. - OBSERVED AND PREDICTED VARIATION OF 
QUENCHING DISTANCE WITH INITIAL TEMPERATURE 


Propane in 
air^ percent 
by volume 

m 

(observed^®' 

d«T-“) 

(predicted. 

Dig 

(predicted. 

3.00 

0.90 

0.54 

1.06 

3.25 

.75 

.49 

1.00 

3.50 

.68 

.38 

.90 

4.03 

.60 

.34 

.85 


^ata from ref. 9. 


Square of quenching distance, d , cm' 



^2 [i-10 ^ p^(ei)/pj 

Figure 1. - Typical correlation of quenching- distance data with consideration of inconplete 
destruction of active particles at vail. 
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(a) Experimental; data from reference 9. 

Figure 2. - Observed and predicted effects of initial temperature 
on quenching distance. 
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Initial temperature, °K 
(b) Predicted) Doo T^. 


Figure 2. - Continued. Observed and predicted effects of Initial 
temperature on quenching distance. 
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Initial temperature, °K 


(c) Predicted; D 

Figure 2. - Concluded. Observed and predicted effects of initial 
temperature on q.uenching distance. 



1 

.2 .4 .6 .8 1 2 4 6 eiD 20 40 60 100 200X10-12 



.1 .2 .4 .6 .8 1 2 4 6 8 10 20 40 60 100 200XL0" 

(a) D OD T^. 

Figure 3. - Correlation of quenching- distance data by means of equation (9). 
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Figure 3. - Concluded 


Correlation of quenching- distance data hy means of equation (9) 
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